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Abstract 

The classical nonlinear laser-plasma interaction theory is corrected. Given the effects of vac- 
uum polarization (induced by extreme laser) as nonlinear media response, one-dimensional wave 
equations of a monochromatic laser field are derived from the Heisenberg-Euler Lagrangian density 
and a derivative correction with the first order quantum electrodynamic (QED) effects. A more 
suitable model to formulate the interactions of extreme laser and high-energy-density plasma is 
developed. In the results, the enhanced effect of vacuum polarization will be discussed and shown. 
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I. INTRODUCTION 



With the growing advancement of laser technology, the power of the next-generation 
optical laser systems will reach about 10 ~ 100PW levels, and after focusing to lfim 2 
the intensity of the laser beam is expected to exceed 10 25 W / /cm 2 [l]. Meanwhile, with the 
electron densities of the plasma approaching unprecedented 10 26 cm -3 achieved under the 
present experimental conditions, X-rays with high energies of E > lOkeV are required to 
penetrate through the plasma^. Taken into account the effects of quantum vacuum in the 
interactions of extreme laser with high-energy-density plasma, the classic theory should be 
modified accordingly. 

As fundamental concepts in quantum field theory, the effects of quantum vacuum were 
neglected for many years but are now beginning to be hot issues for the improvement of 
experiment techniques. The effective Lagrangian density for quantum vacuum corrections 
in electromagnetic field was presented by Heisenberg and Euler as early as 1936 [3j, and 
improved by Schwinger in 1951 [4J. Significantly observable phenomena of quantum vacuum 
such as vacuum breakdown, electron-positron pairs created from vacuum, need to impose a 
super-strong background field. It is the Schwinger critical field E$ c = mlc 3 /eh ~ 10 ie V/ cm, 
corresponding to the laser intensity I cr = E 2 r /8ir = 4.6 x 10 29 W/cm 2 . However, for the 
present experimental conditions it is still very difficult to achieve the critical field-strength, 
and the realization of quantum vacuum lacks experimental validations. Recent years, many 
researchers tried to put forward reasonable models and observable effects which focus on the 
effects of vacuum polarization induced by intense laser pulses {5, (J. The vacuum-polarization 
effects, requiring laser intensity much lower than the critical value I cr , are expected to 
be observed in the near future, and the predictions of QED theory on quantum-vacuum 
phenomena will be verified. In all the models about the effects of vacuum-polarization, the 
interaction of intense laser and plasma is an effective issue. Such research exhibiting the 
effects of QED vacuum via plasma characteristics has indeed made some progress in recent 
years [7-9]. 

When we investigate the processes of intense laser-plasma interaction, an important pa- 
rameter of the laser intensity ao is particularly noteworthy. Defined as the peak amplitude 
of the laser field, normalized vector potential a = eA/m e c 2 , the parameter ao and the peak 
laser light intensity Iq possess the relationship of a 2 , ~ 7.3 x 10~ 19 {\({im)] 2 I (W/cm 2 ). When 
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do <C 1, a linear theory is only required for such interaction lOj, 111). Once in the regime of 
ao> 1, the electronic behaviors are relativistic and the intense laser-plasma interaction must 
be described by a nonlinear theory. A valid one-dimensional analytical model was presented 



in the 1990s by Esarey et al.[12| 
simulation were also obtained 



n the case of three-dimensional, some results via numerical 
14| and subsequent issues related to nonlinear laser-plasma 
interaction, e.g., wave breaking, electron acceleration, multiple laser beam and a series of 
problems were discussed. If the laser light intensity is further raised to make Oq 3> 1, new 
phenomena for the interaction of laser (into the regime of extreme light sources) with plasma 
cannot be explained by the classical nonlinear theory. So it is necessary to complement and 
extend earlier works. 

The effects of quantum vacuum induced by the laser with extreme intensity and high 
frequency can not be neglected. Proper corrections due to the vacuum-polarization effects 
are considered in the coupling equations for describing the interactions of extreme laser and 
high-energy-density plasma. In this paper, the vacuum-polarization effects, similar to the 
response in media, are introduced via Heisenberg-Euler (H-E) correction term and derivative 
correction term in the Lagrangian density of electromagnetic fields. A ID nonlinear cou- 
pling model and the corresponding results are obtained for a more accurate description. The 
characteristics of these two kinds of corrections are compared, and the results are discussed 
in detail with the parameters of the most likely areas, such as inertial confinement fusion 
(ICF) as well as the compact object as dwarf star, for the valid observation in high-energy- 
density plasma physics. Compared with the classical theory of ID nonlinear laser-plasma 
interaction 12] , our model with the contribution of vacuum-polarization effects is more ef- 
fective approach to describe the practical situation of extreme light and high-energy- density 
plasma. 

The paper is organized as follows. In Sec. II, introducing the correction terms of vacuum 
polarization into the equations of electromagnetic field, a ID wave equations of a monochro- 
matic paraxial laser field is derived. In Sec. Ill, a nonlinear coupling model to formulate 
the intense laser-plasma interaction under the influence of vacuum-polarization effects is 
obtained. Finally in Sec. IV, the influence of the vacuum-polarization corrections is in- 
vestigated, and the effective way to enhance the effects in high-energy-density physics is 
explored. 
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II. BASIC EQUATIONS OF LASER FIELD WITH QED VACUUM EFFECTS 



A set of wave equations for laser field has been derived to formulate the extreme light 
behavior. In such wave equations, two kinds of QED vacuum effects have been taken into ac- 
count, (i) The high intensity effect of vacuum polarization is formulated by H-E Lagrangian 
density when the laser intensity is extreme but still much lower than Schwinger critical inten- 
sity, and (ii) the effect of rapid field variations can be illustrated by the derivative correction 
Lagrangian density when the laser frequency is ultra-high but still much lower than Comp- 
ton frequency u e = mc 2 /h. For all these cases, the Lagrangian density of electromagnetic 
fields can be represented as follows [3], 



C = C + SC HE + SC D - -fA 

c 
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where C is the classical Lagrangian density, 5 Che is H-E correction (composed of two 
Lorentz invariants, F^F^ = -2(E 2 — B 2 ) = -4F and F^F^ = -4E • B = -4Q), and 
8Cd is derivative correction. F^ v = d^A v — d v A^ is the four-dimensional electromagnetic 
tensor, F^ u = e^ uXT F\ T /2, □ = d a d a , and is four-dimensional current density, where A 11 
is contravariant electromagnetic potential. The parameter k = a 2 h 3 /457rm e 4 c 5 gives the 
H-E coupling, the parameter a = 4ah 2 /157rm e 2 c 5 gives the derivative effects of vacuum 
polarization, and a = e 2 /he ~ 1/137 is the fine structure constant. 

From the Euler-Lagrange equation, dC/dA^ — d u [dC/d(d u A^)] = 0, the field equations 
with vacuum-polarization corrections can be obtained as 



I + ctU)0„F"" ■ i;<)„ ( l77-'"'' + 1QF" 1 ' ) = ^f, 



(2) 



i.e., 
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(3) 
(4) 



where 5n = n — no is the density perturbation of electron, no is the ambient plasma density, 
and v is the electron fluid velocity. Due to the high intensity effect of vacuum polarization, 



the effective vacuum "sources" (vacuum "current" p vac and vacuum "charge" j vac ) appearing 
in Eqs.flH]) and (J1J can be written as 

Pvac = V- -^(4^E + 7^B), (5) 

47T 

}vac = <%-^(4FE + 7GB) -cVx -£-(4FB - 7£E). (6) 

47T 47T 

Assuming light beams to be monochromatic and paraxial, and neglecting self-focusing 
and diffraction effects in the beam transmission, the propagation of laser pulses can be char- 
acterized by ID model, i.e., V —> e z d/dz = e z d z (for light propagating along z axis). Using 
the Coulomb gauge, V • A = 0, and considering the strong laser pulses with large amplitude, 
the longitudinal component (the one perpendicular to z axis) of the vector potential is much 
larger than the transverse one, |A^| <C |Aj_|, and then A « A^ = v4e_|_. The effective 
polarization P mc and magnetization M mc induced by vacuum "sources" are written simply 
as 

Pvac = -—Fd t Ae ± , (7) 

TIC 

M mc = — Fd z A(e z x ej_). (8) 

7T 

In the ID model, the electromagnetic potentials of the laser beam are demonstrated as 
A(z, t) = A (z, t) e** and $(z, t) = %(z, t) e m , where A (A , O ) and * = -ut+kz+ip(z, t) 
are the amplitude and the phase of the laser potential, respectively. The perturbations of the 
amplitude 6Ao(6A , <50 o ) and the phase Sipo induced by the vacuum-polarization effects are 
much small- varying, i.e., \8Ao\ <C |^4o| and \5ip\ <C |^|. So the first-order differential terms of 
the amplitude and the phase will only be kept when we calculate the corrections of vacuum 
"source", and the approximate relations are derived as diAdjA « (d{djA)A, di<ftdj<f) ~ 
(didj4>)4>, and diAdj(f) ps {didjA)4> ps (didj(f))A ps di<fidjA, where i,j = t,z. 

Since the H-E coupling coefficient k is much small (k ~ 10~ 32 cm ■ s 2 /g) and the laser 
field is much lower than Es c , the intensity contribution of vacuum-polarization correc- 
tions (vacuum "source") can be written as a perturbation term. Neglecting the higher- 
order terms {0(k)) in the equations, the vacuum "source" are simplified as follows, 
Pvac ~ —6k • 47re (5n$ — n(3±A) e5n and j vac ~ 6k • Ane (5n$ — n(3±A) en(v z e z + v j_e_i_), 
where 5n ~ <9 Z 2 $, n/3j_ ~ (d z 2 — d t 2 /c 2 )A. Similarly, considering the coupling coefficient a 
of the derivative correction is small (about 10~ 24 cm 2 ) and the laser frequency is much lower 
than u e , the first-order term of a is only kept in the equations. Finally from Eqs.(j3]) and 



(jl]), the wave equations of a monochromatic paraxial light beam with vacuum-polarization 
corrections can be represented as 

(V - ^ d t 2 ) a = V^-L - 6«:(47re) 2 (5n0 - n^ ± a)n/3 ± - k 2 a (d 2 - ^d t 2 ^j (j^P±^ , 
-d t d z( f> = -k 2 —(3 z + 6K(4ne) 2 (6n4> - n(3 ± a)nf3 z + k 2 a id 2 - \dA (~M , 

d 2 <f) = k 2 — - 6/<47re) 2 (cW> - n(5 L a)5n - k 2 a (d 2 - ( —) , 

' n p \ c 2 J \n J 



(9) 



where = (e/m e c 2 ) $ and a = (e/m e c 2 ) A are the normalized electromagnetic potentials, 

k 2 = lu 2 /c 2 = 4iTe 2 n /m e c 2 , 0± = v_l/c, and (3 Z = v z /c. 



III. THE RESPONSE OF VACUUM POLARIZATION IN THE NONLINEAR IN- 
TENSE LASER-PLASMA INTERACTIONS 

In the regime of extreme lasers, the processes of intense laser-plasma interactions are 
something different from the classical theories. Especially when the intensity of laser beams 
is approximated to the Schwinger critical intensity, the QED vacuum effects will play a 
significant role. In this section, considering the proper corrections of vacuum polarization, a 
ID nonlinear coupling model is found to describe the self-consistent interactions of intense 
laser with plasma . 

To represent the plasma wave (laser wakefield) generated by this interaction, the fluid 



plasma model is generally used 



15| . The perturbation of the plasma electrons (plasma wave) 



generated by the laser fields obeys the following hydrodynamic equations, 

d t n + V • (nv) = 0, 

d t p + (v ■ V)p = -e (E + ^ x b) , 

where p = mv7 is the momentum of electron, 7 = (1 — (3 2 ) is the relativistic factor and 
P = v/c. 

Taking account of the vacuum-polarization effects, a driving laser pulse with extreme 
intensity can be formulated by Eq.Q. Assuming that the laser wakefield is under the condi- 
tions of paraxial approximation, the velocity of the electrons in wakefield can be decomposed 
into transverse and longitudinal, (3e = (3±e± + (3 z e z , in ID model. So the continuity equation 



6 



IS 

^d t n + d z (n/3 z ) = 0, (10) 

and the momentum equation of electronic fluid in laser wakefield is divided into two com- 
ponents (transverse and longitudinal), 

\d t {lPi) + PMlPd = - c d t a + (3 z (d z a), (11) 

and 

\dt(rtPz) + MzWz) = d z <p - 0x(d z a). (12) 
The transverse component, perpendicular to the direction of wake- field propagation, 
represents the rapid quiver motion of plasma electrons under the disturbances of intense 
laser pulse. From Eq.( lTT|) . normalized vector potential a = 7/3j_ and relativistic factor 
7 2 = 1 — P 2 — Pi = (1 + a 2 ) I (1 — P z 2 ) can be easily obtained. 

The longitudinal component mainly formulates the drift motion of electrons, the accel- 
eration or deceleration motion in the plasma wakefield, induced by the intense laser along 
its propagation direction. The vacuum-polarization effects fluctuate the plasma wave by 
influencing the driving laser directly. 

When the intensity of driving laser beam exceeds a certain value (about 10 18 W/cm 2 for 
Ifim wave length), i.e., a ~ 1, the electron motion in the laser wakefield is relativistic. 
Further, when a ^> 1, the longitudinal drift velocity of electrons (approaching the speed 
of light) is much larger than the transverse quiver velocity, namely, (3± <C /3 Z ~ 1. From 
Eq. ffl2l . the fluid equations can be obtained in the form 

\<kP z = --^ (d z + - c P z dA a 2 + 1 (1 - P 2 ) d z <j>, (13) 

-da = P z d z <P + ^-d z a 2 , (14) 
c 2c7 

where d t — d t + cP z d z . It is convenient to perform a transformation from the laboratory 
frame to the speed of light frame, £ = z — ct, t = t. The potential wave Eq.([9]) of the intense 
laser beam can be rewritten as 

n 
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dp 2 (f) = & 2 — — 6ft(47re) 2 ( 
s n V 



5n(f) — —a 2 ] 5n — k 2 a \-dtd T — \d 2 J — , 

7 J p \c & ) n 



and 



[7(1 - Pz) 



-d T n = d ( [n(l - P z )\ 
-0] = --3 T ( 7 &). 



Assuming u; ^> w„, we can use the quasistatic approximation, 



d T / nd£' 



(16) 
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< 



no, 



<C 1, and obtain the plasma wave hydrodynamic equations, 
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(18) 
(19) 
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Representing the intense laser-plasma interactions with the contribution of the vacuum- 
polarization effects, the coupling Eq.(fl5|) can be simplified as 



dfr = -fv 



1 + a 2 

:i+<t>) 



(21) 
(22) 



and 



1 - A HE - A D 

(a 2 + 2 )(2 + 0) 



1 - A 6 n 



(23) 



where 77 is the correction factor, A 3 = (4/15)a 2 A c 3 is the coupling coefficient of vacuum- 
polarization , and X c = h/m e c is the reduced Compton wavelength. Eqs.f l2Tj) and fl22|) 
are the expressions of the ID nonlinear coupling fields containing the contribution of the 
vacuum-polarization effects. Eq. (|2ip mainly formulates the propagation characteristics of 
extreme laser through the plasma, and Eq. (1221) is about the evolution of plasma wave induced 
by the extreme laser beam. The first term on the right of Eq. ([2"B"j) is the classical one, the 
second term Ahe is attributed to the H-E correction and the third term to the derivative 
correction. Without considering the contribution of vacuum-polarization effects, i.e., A HE ~ 
A D — > 0, fj ~ 1, the Eqs. (l2"T|) and fT2~T]) reduced to the nonlinear classical equations founded 
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by Esarsy et al. 12j, |l3j. This coupling model with the vacuum-polarization effects can 
help us to investigate the transmission characteristics of laser-plasma interaction under the 
extreme field laser. 



IV. CONCLUSIONS AND DISCUSSION 

In this work, a ID nonlinear model of laser-plasma interaction containing vacuum- 
polarization corrections has been studied. This nonlinear model of the coupling field de- 
scribes the impact of extreme laser on plasmas. We will focus the attention on how to 
significantly enhance the effects of vacuum polarization in the coupling processes as follows. 

Taking a 3> 1 into account, the H-E correction rate, Ahe — 
A 3 no(a 2 + 2 )(2 + 0)/(l + 0) 2 , can be simplified as A 3 n a 2 . According to the afore- 
mentioned relationship of a and /, it is easy to obtain the relation as n^a 2 ~ 816(a;p/a;) 2 /. 
Assuming u p /u> ~ 0.1 and / ~ 10 25 W/cm 2 , the H-E correction rate Ahe is closed to 1CP 11 . 

The derivative correction rate, A D = A 3 n (47r//c 2 )(2c _1 <9,c<9 T — c~ 2 d T 2 ), is reduced to 
A 3 47rx(w/u;p) 2 no under the quasi-plane-wave approximation, where x = 1 — k 2 c 2 /u> 2 is the 
magnetic susceptibility of plasma determined by specific systems. Assuming that u p /u ~ 
0.01, x ~ lj an d n Q ~ 10 25 cm 3 , the derivative correction rate is Ap ~ 10~ 6 . 

The H-E correction of vacuum-polarization is decided by the laser intensity directly. It is 
calculated that the H-E correction rate A^e is no more than 10~ 10 when laser beams with 
/ ~ 10 25 W/cm 2 travel in an underdense plasma. In a dense plasma, the H-E correction 
is weakened but more complicated and worthy for further studying, whereas the derivative 
effect becomes obvious. Especially by using driving laser beams with ultra-high frequency, 
the observation of the vacuum-polarization effects is expected to be implemented. Figs. [T] 
and [2] show the relationship between the rates of the vacuum-polarization corrections and 
the important parameters in the laser-plasma interaction. 

In future ICF experiments, which are laser-plasma interactions with the highest plasma 
density in the experimental environment, electron density no approaching 10 -26 cm -3 is 
expected to become possible. For this extremely dense states, the derivative correction 
rate is An ~ 10~ 5 . In astrophysics, with n ~ 10 -30 cm -3 in the compact objects such as 



dwarfs 



16| , the magnitude of Ap can reach 10 1 further. 
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FIG. 1: The magnitude of derivative correction rate in different plasma density uq for \ ~ 1) 
/ = 10- 25 cm" 3 and uj p /oj = 0.01. 
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FIG. 2: The magnitude of vacuum-polarization correction rates Ahe and Ad in different frequency 
factor uj p /oj for x ~ 1, I = 10~ 25 W/cm 2 and n = 10 25 cm~ 3 . 
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